Statistical properties of phases and delay times of the one-dimensional Anderson 

model with one open channel 



> 

o 



CZ2 



a 



X3 

o 



> 

oo 

o 



o 



X 



A. Ossipov, Tsampikos Kottos, and T. Geisel 

Max-Planck-Institut fiir Stromungsforschung und Fakultat Physik der Universitat Gottingen, 

BunsenstraBe 10, D-37073 Gottingen, Germany 

We study the distribution of phases and of Wigner delay times for a one- dimensional Anderson 
model with one open channel. Our approach, based on classical Hamiltonian maps, allows us an 
analytical treatment. We find that the distribution of phases depends drastically on the parameter 
aA ~ a /sink where a^ is the variance of the disorder distribution and k the wavevector. It under- 
goes a transition from uniformity to singular behaviour as a a increases. The distribution of delay 
times shows universal power law tails 1/r^, while the short time behaviour is a a- dependent. 
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Quantum mechanical scattering by chaotic or disor- 
dered systems, has been a subject of a rather intensive 
research activity during the last ten years. This inter- 
est was motivated by various areas of physics, ranging 
from nuclear M^ atomic and molecular |^ physics, 
to mesoscopics y] and classical wave scattering [^ . The 
most fundamental object characterizing the process of 
quantum scattering is the unitary 5'-matrix relating the 
amplitudes of incoming waves to the amplitudes of out- 
going waves. At present, there are two complementary 
theoretical tools employed to calculate statistical prop- 
erties of the 5— matrix, namely the semiclassical and the 
stochastic approach. The starting point of the first is a 
representation of the S— matrix elements in terms of a 
sum of classical orbits MM while the later exploits the 
similarity with ensembles of Random Matrices [[7|. Thus, 
for chaotic systems, many results are known. Namely, the 
distribution of the 5— matrix was found to be described 
by the so-called Poisson Kernel |g|,0 . In the same frame- 
work, the Wigner delay time statistics has been studied 
intensively. This quantity captures the time- dependent 
aspects of quantum scattering. It is related to the time 
spent in the interaction region by a wavepacket of energy 
peaked at E. The delay time is not self-averaging and 
one must have its full probability distribution. For the 
one-channel case this was found in H while the case with 
generally M open channels was presented in g . 

In spite of the wealth of new results on chaotic scatter- 
ing, very little is known for the scattering from systems 
being in the localized regime P], even for the simplest 
case of one dimensional (ID) disordered systems where 
Anderson localization dominates [10[ . For the latter case 
and for only one channel, there are some analytical re- 
sults about the distribution of phases of the S— matrix 
and of delay times. They were first derived with the 
help of the invariant embedding method pd| ] and more 
recently with the use of exponential functionals of Brow- 
nian motion [Q. In all cases, however, these results are 
restricted to continuous systems and to weak disorder. 

In this letter, we surmount these limitations by using a 



different approach, which transforms the initial ID tight- 
binding Anderson model into a classical linear oscillator 
with a parametric perturbation given in the form of peri- 
odic (5- kicks ]l^,|ljl . The amplitudes of these kicks are de- 
fined by the ratio An = Vn/ sink, where Vn is the on-site 
random potential with variance cr^ and k is the wavevec- 
tor. Based on this approach we derive simple iteration 
relations for the phases of the S'-matrix and the delay 
times which allow us to analyze their distributions. We 
find that the kick strength a a = a j sink- and not o itself 
pG-O]- is a new significant parameter that controls the 
form of these distributions . Wc challenge - even for the 
weak disorder regime ct <C 1 - the one-parameter scaling 
based on the so-called random phase hypothesis |ll5| , and 
conclude, in agreement with pM, for a new criterion for 
phase randomization which is cTyi <C 1. In this limit, we 
derive the distribution of delay times analytically. It is 
found to be insensitive to the specific type of the random 
potential, but may depend on its variance. On the other 
hand, in the opposite limit cr^ ^ 1, it becomes disorder 
dependent. The tail however follows the same power law 
1/t^, as in the cr^ "C 1 limit. 

We consider a ID disordered sample of length L with 
one semi-infinite perfect lead attached on the left side. 
The system is described by the tight-binding equation: 



V'n-t-i + V'«-i = (-E- K)V'«; E^2cosk 



(1) 



where ipn is the wavefunction amplitude at site n. We as- 
sume that for < n < L, Vn is random delta-correlated 
given by a distribution Vv with mean zero and variance 
(T^ . For n < 0, V^ = and we impose Dirichlet boundary 
conditions at the edge V'i+i = 0- Therefore, for n < 0, 
scattering states of the form ijjn = e**^" -I- re"'*^" rep- 
resent the superposition of an incoming and a refiected 
plane wave. Since there is only backscattering, the re- 
flection coefficient r{E) = e**^-^) is of unit modulus and 
the total information about the scattering is contained 
in the phase $(-E). The Wigner delay time is given by 
t{E) = ^^ ' . Our aim is to find the probability distri- 
bution of the phases 7'$(<I') and of delay times "Pr (''')• 



Eqn. (Wj can be written equivalently in a form 



Consequently we have the following relation for the total 



transfer matrix P = Jln=o -^^" 



Pll Pl2 



i^L ) V P21 P22 I \ e~^^ + re*'= 



1 + r 



Solving Eqn. (0) for r we get r 
where </> = $/2 is now given by 



Pll+Pl2C~"' 
Pll+Pl20''= 



tan(j) = 



Pll + Pl2C0sk 

Pi2sink 



(3) 



c,2i4 



(4) 



As was indicated in UM, one can rewrite (0) in the form 
of a two-dimensional Hamiltonian map Q„ 



x„+i \ n I ^" \ n I cos/c — An-sink sink 
Pn+i J ^^ \ Pn J ' " \ —Ancosk — sink cask 



Here, (x„,p„) play the role of a position and momen- 
tum of a parametric linear oscillator subjected to periodic 
kicks of strength An = j^ and period T = 1. Between 
two successive kicks, there is a free rotation in the phase 
space which is determined by the eigenenergy E of our 
initial equation (|l|). In such a representation, amplitudes 
■i/'n of a specific eigenstate correspond to positions of the 
oscillator at times t„ = n. Qn is related to the transfer 
matrix Af„ through a similarity transformation p3[ : 

Qn^RMnR-'; R={,,J/,,^k-l/s^nk)■ («) 

As a result, the total transfer matrix P is related to the 
map F = nri=o *5" through the similarity transforma- 
tion P = R^^FR. Using this, together with (^ we write 
(j) in terms of our Hamiltonian map F as: 



tancj) 



111 

-F, 



12 



We can give a geometrical interpretation for (0). Con- 
sider the time evolution of the vector v{t = 0) = (0, 1)"^ 
under the inverse map F~^. For time tn = L we have 



v{t ^ L) = P- 



1^0 



-F12 
Fii 



(8) 



Then ip is exactly the angle between the vector v{t — L) 
and the x-axis. It is convenient therefore to pass to polar 
coordinates (r„,0„), using the transformation x = rcosO 
and p = rsinO. Then (ph is written in the form [ p3[ : 

^l+i = rlDn -,01 = {1 + Alcos^0n - Ansin2en) 
cos6n+i = Dn^[cos{9n + k) + AnCosOnSink] 
sinOn+i = P>n'^[sin{6n + k) — AnCOsOnCOsk]. (9) 



The relation between (^„+i and 0„ is found after in- 
version of (0) and has the form IIitI: 



tan{(t>n+i) = tan{4>n - fc) + A„ 



(10) 



From (lOp, we get the following equation for the sta- 
tionary (n — > cx)) distribution of phases 7^^ ((/)): 

V^iatany) ^ j" ^ V^{atan{y ~ s) + k) 



i + r 



1 + (2/ - sy 



In the cTyi <C 1 limit, we can approximate Va{s) by a 
i5— function. Then from { \\/\\ ) we have for the 'P($) 



V^i<^) = V^i<^ + k) 



(12) 



which for k equals to irrational multiples of tt, leads to a 
uniform distribution 



7'*($) 



1 

2^' 



k^Q 



(13) 



(5) 

in agreement with previous numerical results |10[| . We 
understand (13) in the following way: For cr^ <C 1, the 
particle travels long distances inside the sample and un- 
dergoes many scattering events which leads to random- 
ization of the phase. We should stress that for k equals 
to rational values of tt, the distribution is non-uniform. 
For example, for k — 7r/2 we obtain |14 17| , p^ 



P*(*) = 



1 



2^-)V^ 



cos(2$; 



k = 



(14) 



where /C is the complete elliptic integral of the first kind. 
A similar anomaly at the band center governs the dis- 
tribution of the angles 6 of the direct map (H) [Q . The 
above expressions ( |l3| , |l4| ) describe well the corresponding 
numerical results presented in Figs. la,b. 

In the a A 3> 1 limit, we distinguish between two parts 
of the spectrum. Namely, k near the band edges (i.e. 
fc ~ 0, tt) and far away from them. For the latter case 
(7) and for Vn distributed uniformly between [— Y't] '-^^^ 
can derive an analytical expression for 7^$ by rewriting 
( |10| ) in terms of the probability distributions Vtan{4>)- We 
get (in first order in 1/<ta) the uniform distribution 118] 



Vt. 



.Wiy) = (\/l2^A)"'e (VSaA -\y- cotk\^ . (15) 



where Q{x) is the Heavyside function. The distribution 
of phases ■P$($) can be derived easily from (|l5|): 



7'*($) = 



0.5 Q (vSaA — \tan{^) — cotk 



V3< 



(^A 



C0s2(f) 



(16) 



In contrast with the cr^i ^ 1 limit (n3[) , here the distribu- 
tion of phases (Hq) is highly non-uniform. By increasing 
the disorder strength u, two peaks appear in the neigh- 
bourhood of $ = TT (see Fig. Ic) while a gap is created 



between them. As cta increases further, the two peaks 
move closer to one another. Now the scattering is so 
strong that most of the particles are reflected back from 
near the surface, thus being scattered only from a few 
sites. In the limit cr^ ^ oo, the distribution becomes a S 
function centered at <I> = tt. ( p^ ) is in very good agree- 
ment with our numerical data presented in Fig. Ic. The 
behavior near $ = tt is more subtle than the one given 
above; nevertheless, (H) gives the correct scale on which 
the distribution vanishes near ^ — n. Distribution dlq) 
agrees quite nicely also for other disorder potentials pSl . 
Near the band edge (fc w 0, tt), our detail numerical 
analysis showed that 7'$($) is again highly non- uniform. 
Namely the distribution becomes narrower and centered 
(although with a slight asymmetry) at ^ — tt . We no- 
tice that such a choice of the parameters, can be realized 
even for weak disorder a <C 1. In Fig. Id we present a 
representative case corresponding to fc = 10^"^ -^/tF and 
a = 0.0577. 
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FIG. 1. Distribution of the phases P($) for various disordered 
strengths and energies. The on-site potential is uniformly dis- 
tributed between [--j; -y]- (a) V = 0.2 and A; = 0F. (h) V = 0.2 
but now for k = tt/2 (band center); (c) V = 10 and k = ^/tt. Our 
numerical data (histogram) are in perfect agreement with the an- 
alytical predictions (dashed lines), (d) V = 0.2 and k = 10^''v''i"- 

Let us now turn to the analysis of delay times. From 
( [lO| ) we get the following iteration relation for t„ 



Tn+l 



G 



1 



sink 
1 + Ansin (2((/)„ 



A„ 



cotk 



1 -I- {tan{<pri 



k)) 



Alco.'Pl 



k) + An) sink 
.-k) (17) 



which proves to be very convenient for numerical calcu- 
lations since it anticipates the numerical differentiation 
which is a rather unstable operation. 

For cr^ ^ 1 and k equals to irrational multiples of tt, 
we obtain an analytical expression for Vr{T). To this end, 
we first derive from (HG) and (O) a stochastic equation 
for the phases 4> and rescale delay times f = ct^t (up to 
a'^j^, respectively p^: 



—- ~ — fc — aj^sin{(j) — k)cos {(p — k) + cos {(j) — k)A 



dt 
df 

'di 



— ~ — cr^ (f(cos {(p — k) — sin 2{<j) — fc)) — sink) 

-fsin{(l)-k)A. (18) 



Using the van Kampen lemma, from (h8[) we obtain the 
Fokker-Planck equation for the joint probability distribu- 
tion of (j) and f. Using the fact that (p follows the uniform 
distribution (see Eqn. dl3)) and assuming that the vari- 
ables (j) and T are statistically independent, we obtain the 
Fokker-Planck equation for Vf (f) 



dt 



o\ 






■isink)Vf{f,t)) 






(19) 



The resulting stationary distribution is obtained by set- 
ting ^^ = and has the form 



V{t) 



h 



VT^ 



^expi-loo/vr); cta < 1 



(20) 



where loc = 2(4 — E'^)/a'^ is the localization length and 
V — \dE/dk\ is the group velocity. We note that ( poj ) is 
independent of the nature of the disorder; it only depends 
on its second moment through the localization length. 
( P0| ) takes its maximum value at Tmax = 0-5loo/v, in- 
dicating that the most probable trajectory that an elec- 
tron travels (forth and back) before it scatters outside the 
sample is the mean free path l]\i — /oo/4. Our numerical 
results (see Fig. 2a) are in perfect agreement with (|2C|). 
An expression similar to ( pO| ) was obtained in |ll| and 
more recently in |]12[ using totally different approaches. 
It is important to stress here that the latter analytical 
results refer to continuous disordered models. 

Finally we discuss the distribution of delay times V{t) 
for (T^ ^ 1. In this limit, we were not able to derive any 
analytical expression. Our iteration relation (p^ ) how- 
ever, has proven very efficient for numerical investiga- 
tions. In Fig. 2b we show the distribution of the delay 
times for a uniform and a Gaussian on-site potential dis- 
tribution with the same variance a^ — 10. It is clear that 
the short time distribution differs considerably in the two 
cases and also with respect to the theoretical prediction 
(p0|). On the other hand, the distribution of large delay 
times, reflecting the times for which the wave penetrates 
deeper into the sample shows the same 1/t^ behavior in- 
dependent of the form of the disorder potential. It is this 
part of the reflected wave, and not the prompt part that 
is expected to behave in a universal way p2[ . The asymp- 
totics of the distribution for large r is presented in the 
inset of Fig. 2b, where we plot the integrated distribution 
7(7) = Jq P{Y)dY of the inverse delay time 7 = 1/r in 



log— log scale. To this end, we calculated 10^ delay times 
using the iteration relation (|l^. In both cases (Gaussian 
and uniform distribution) presented in the inset of Fig. 2b 
we collected at least 10** delay times that were larger than 
T > 150. Our numerical data clearly show that 1(7) ^ 7 
for 7 <C 1. Thus 7'(t) ^ l/r^ in agreement with the 
above mentioned universal behavior |1Q,^31 . 




10000 



to a logarithmic divergence of the average value of r indi- 
cating the possibility of the particle traversing the infinite 
sample before being totally reflected. As was indicated 
in ||ll|,|ljl this is another manifestation of the fact that 
in the localized regime the conductance shows lognormal 
distribution due to the presence of Azbel resonances. 

We acknowledge useful discussions with F. Izrailev and 
I. Guarneri. (T.K) thanks U. Smilansky for initiating his 
interest in quantum scattering. 
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FIG. 2. (a) Distribution of the delay times "Pir) for on-site po- 
tential, uniformly distributed between [—0.1; 0.1] and wavenumber 
k = y/n. The dashed line corresponds to (EOl). (b) Distribution 
of the delay times VIt) for uniform and Gaussian Vv- In both 
cases the variance is a^ = 10 and the wavenumber is fc = ^/n. The 
dashed line corresponds to (M). The inset shows the integrated 
distributions of inverse delay time 7 = l/r. The dashed line has 
slope 1 and is drawn to guide the eye. 

To summarize, we analyzed the distributions of phases 
7^($) and delay times V{t) for a ID Anderson model with 
one channel in the localized regime. For V{^) we found 
that it undergoes a transition from uniform {a a < 1) to 
non-uniform behaviour (a a > 1), where a a — a /sink. 
The latter limit can be achieved either by decreasing the 
disorder strength a, or by taking fc — > 0, tt. This has fur- 
ther implications for the so-called one-parameter scaling 
hypothesis. The scaling theory of disordered systems |1^] 
assumes that V{^) is always uniform in the weak disorder 
limit tr <?C 1 leading to the one parameter scaling. Viola- 
tion of the latter is always associated with an increase of 
a. Our analytical and numerical results showed however 
that a strong violation of the phase uniformity occurs 
even in the case of weak disorder ct <C 1 , in an apparent 
contradiction with the existing picture. In this limit, the 
states near the band edge fc « 0, tt never obey the single- 
scaling hypothesis since cr^ :^ 1. Thus the spectrum 
of the system is divided into two groups with different 
scaling properties, which coexist at the same strength of 
disorder |16|. For V{t) we found that only the short time 
behaviour is affected as we are increasing a a- For t —^ 00, 
V{t) ~ 1/r^ independent of the value of ct^. This leads [18 
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